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Abstract 

A two dimensional flat phononic crystal (PC) lens for focusing off-plane shear waves is proposed. 
The lens consists of a triangular lattice hole-array, embedded in solid matrix. Self-collimation 
effect is employed to guide the shear waves propagating through the lens along specific directions. 
Dirichlet-to-Neumann Maps (DtN) method is employed to calculate the band structure of the 
PC, which can avoid the problems of bad convergence and fake bands automatically in the void- 
solid PC structure. When the lens is illuminated by the off-plane shear waves emanating from a 
point source, a subwavelength image appears in the far-field zone. The imaging characteristics are 
investigated by calculating the displacement fields explicitly using the multiple-scattering method, 
and the results are in good agreement with the ray-trace predictions. Our results may provide 
insights for designing new phononic devices. 



PACS numbers: 41.20.Jb, 43.20.E1 



Phononic crystals (PC) are artificial periodic elastic media made for controlling the prop- 
agation of acoustic/elastic waves. The propagating modes of the elastic waves in PC are 
Bloch waves, and their dispersion relations form the phononic band structures, containing 
passbands and bandgaps. Various structures for reflecting, trapping, and guiding acous- 
tic/elastic waves relying on the existence of PC bandgaps have been proposed and studied 
[IH9]. Besides the bandgaps, the phononic passbands are also useful and have interesting 
properties. When operating in the linear region well below the top of the first band (i.e., 
the long- wavelength region), a PC behaves just like a homogeneous medium, and lens-like 
devices can be constructed from it [TUl E]. Recently, the study of wave propagation behav- 
iors beyond the long- wavelength region has attracted a great deal of interest. It has been 
shown that a slab of PC operating at an appropriate frequency can focus the acoustic waves 
emanating from a point source to form a subwavelength image, overcoming the diffraction 

limit pans!. 

To achieve the slab imaging ability, either negative refraction or anomalous refraction 
is required [T3HTT] . which means the direction of the averaged Pointying vector (S) must 
have a component antiparallel to the Bloch wavevector k or some effective Bloch wavevector 
k e / = k — ko, i,e., k e / ■ (S) < 0. Here ko stands for certain symmetry point in the first 
Brillouin zone or on the zone boundary, such as the M point for the square-lattice structure 
of PC. Recently, it has been pointed out that the PC band structure can also provide other 
wave guiding mechanisms such as self-collimation [TBH2I] and canalization [TTl [22] • Using 
these mechanisms, waves can be guided along certain specific directions inside the slab and 
thus some desired devices such as hyperlens can be realized [23J. 

Up to now most of the studies on the negative refraction or subwavelength imaging 
phenomena adopted the PC that sustain pressure waves, and the discussions on similar 
phenomena in solid elastic media are relatively few [I5HT7] . In this letter, we employ the 
self-collimation effect of PC to design a flat lens for focusing off-plane pure shear wave. The 
lens consists of a triangular lattice hole-array embedded in a solid matrix. When the lens 
is illuminated by the off-plane shear waves emanating from a point source, a subwavelength 
image appears in the far-field zone. The imaging characteristics for a finite lens are inves- 
tigated by calculating the displacement fields using the Multiple Scattering (MS) method 
[TT] |21j . The PC band structure is obtained with the Dirichlet-to-Neumann Maps (DtN) 
method [251 125] . Effective negative refractive index of the PC can be extracted from the 
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equal frequency contour (EFC), and the image points can be predicted according to Snell's 
law. The results are in good agreement with the calculations of field distributions. 

It is known that in calculating the band structure of the fluid-solid or void-solid PC, the 
traditional plane wave expansion (PWE) method encounters the problems of bad conver- 
gence and fake flat bands [271 '28\. These problems can be resolved by replacing the void 
part with a "pseudosolid" or "low impedance material" [291 130] , which is an imagined solid 
material having a very low density and appropriately chosen wave speed (usually very large). 
When these material parameters are properly chosen, the flat band will appear only in the 
uninterested high frequency range, and thus we can safely obtain the desired band structure 
in the moderate frequency range. Although this modified PWE method can indeed avoid 
the flat band problem, there is no definite way of choosing the material parameters of the 
pseudosolid, and usually some experiences are needed. In order to reduce the uncertainty, 
we adopt instead the DtN method to calculate the band structure of the PC. Since this 
method can be applied to the void-solid case directly, the flat band problem disappears 
automatically. Up to now DtN method has only been used in calculating the band struc- 
tures or guided modes of two-dimensional photonic crystal (PhC) systems, however, there 
is no essential difference mathematically between the two-dimensional (2D) electromagnetic 
(EM) and acoustic (AC) systems if the polarizations of the waves are properly chosen. For 
example, the transverse-magnetic (TM) wave propagating in a 2D PhC and the pure shear 
wave propagating in a 2D solid PC are governed by the same kind of scalar wave equation. 
Thus an EM wave problem can be easily translated to an AC wave problem and a one-to- 
one correspondence can be established between these two systems [31]. This consideration 
implies that DtN method can be used in the present case of pure shear wave. A detailed 
description of DtN method will be given later. 

Since in our chosen structure the off-plane shear wave cannot penetrate into the holes, 
the band structure of the PC is determined only by the geometrical factors (the hole size, 
lattice constant, and lattice type) and the material parameters (the mass density and the 
Lame constant for shear wave) of the host medium. In fact, if the lattice type, hole shape 
and the filling fraction of the hole in a unit cell are all given, the band structure of the PC 
becomes universal [33j [31] if the dimensionless frequency ua/27rc t is used as the frequency 
variable. Here c t represents the propagation speed of the shear wave in the host medium 
and a is the lattice constant. If we replace the original host medium by another material, 
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the universal band structure predicts that identical wave phenomena can happen in this 
new PC, usually in a different frequency or for a different PC size. Note that although for 
non-dispersive host media the band structure for a PC consisting of several materials can 
also be expressed in a size-independent universal manner, the material-independence is a 
unique property of the hole-array PC. Moreover, in this case the (shear) wave speed is the 
only relevant material parameter to determine the band structure. 

To be more specific, we choose aluminum as the host medium (p = 2.692g/cm 3 , c t = 
3220m/s). The lattice type is triangular, the lattice constant is a, and the radius of the 
holes is r = 0.4a. An off-plane monochromatic shear wave u(r, t) = zu(r)e~ tuJt propagating 
(in plane) in the host medium satisfies the Helmholtz equation 

V 2 w(r) + k 2 u(r) = 0. (1) 

Here u(r) = u(x,y) represents the amplitude of the displacement, k t = oo/ct is the wave 
number of the shear wave, and the V 2 = d 2 /dx 2 + d 2 /dy 2 is the 2D Laplacian operator 
because u does not depend on z. Since the holes are all circular, displacement field in the 
host medium around a hole (say, the ith hole) can be expressed as a sum of cylindrical waves: 

u(r) = f; C n (j n (k t r) - y^jrY^kr)] e**. (2) 

n=— oo ^ n ^ ' ' 

Here r = |r — is the magnitude of the position vector relative to the center of the ith hole, 
and <p = (pr-r, is its azimuthal angle. To derive this expression, the boundary condition that 
the shear stress vanishes along the hole boundary has been applied. 

Equation ^ says that both u and its first derivative alone any direction are determined 
by the coefficients C n . Therefore, a DtN map between u and its normal derivative du/dn 
along the boundary of a unit cell can be established. Using this relation together with Bloch 
boundary conditions u(r + ax) = e lk ' ai ii(r), u(r + a.2) = e* k ' a2 w(r), the Bloch wave modes 
u(r) and the PC band structure can be determined [251126] . Here ai^ represent the primitive 
lattice vectors, and k is the Bloch wave vector. Equation similar to ^ but expanded as 
sum of Bessel and Hankel functions can also be derived. Using that equation, together with 
the fact that the incident wave around a hole comes from the scattered waves of other holes 
and the waves from the sources, a matrix equation relating the scattering coefficients and 
the source coefficients can be established. Solving this matrix equation then gives us the 
displacement field in the host medium under the illumination of the sources [HI [21] . 
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The band structure and the EFC of the PC calculated using DtN method are shown in 
Fig.l (a) and 1(b). We have also calculated the same band structure using the pseudosolid 
method, and no visible difference between the results of these two methods can be found. 
The accuracy of the band structure is thus confirmed. The two inclined straight lines in Fig. 
1(a) stand for the "sound lines" of the host medium. To focus the shear waves emanating 
from a point source, we choose k t a/2ir = ua/27ic t = 0.6 as the operating (reduced) frequency. 
The EFC of the PC at this frequency looks like a hexagon, which occupies an area a little 
smaller than the corresponding circular EFC of the host medium. The frequency surface 
u;(k) of the PC in the 2nd band around the T point has a negative curvature (negative 
effective phonon mass), which leads to k • \ g < 0, and thus negative refraction can occur. 
Here w g = VkW is the group velocity, pointing in the direction of the greatest rate of increase 
of u and is perpendicular to EFC. It can also be shown that the direction of v s is the same 
as the direction of energy flow (averaged Poynting vector) [32] . When a plane wave of wave 
vector k is incident upon a PC, the refracted wave vector k' can be determined by three 
rules. First, the incident and refracted waves have the same frequency. Thus the tips of k 
and k' are located on the host and PC EFCs, respectively. Second, the wave phase on the 
two sides of the host-PC interface must be the same, which means the tangential component 
of k is equal to that of k'. Third, the normal component of v 9 does not change sign in 
the refraction process, which is a consequence of the energy conservation law V ■ (S) = 
for monchromatic wave. An illustration of this refraction process in terms of the k-space 
description and its corresponding r-space configuration is shown in Fig. 1(b). The thick black 
arrows represent the incident and refracted wave vectors, whereas the pink arrows indicate 
the v 9 (energy flow) directions. Hereafter we always assume TK be the orientation of the 
host-PC interface. The EFC plot in Fig. 1(b) also indicates that if the incidence angle is 
in the interval (— 7r/3, vt/3), wave can penetrate into the PC and the refraction angle of ray 
(energy flow) are almost fixed at ±7r/6. This is caused by the self-collimation effect of the 
hexagonal shaped EFC. On the other hand, if the incident angle is beyond (— 7r/3, 7r/3), 
no refracted wave can be excited, and the incident wave will be totally reflected. These 
phenomena are illustrated in Fig. 2 (a)-(c). Here we use Gaussian beams as incidence waves 
instead of plane waves. The Gaussian source used in these simulations is formed by giving a 
Gaussian strength profile along 20 aligned point sources distributed evenly within 4a. The 
incidence angles in Fig. 2(a), (b) and (c) are vr/18 (10°) and 2vr/9 (40°), and tt/3 (60°), 
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respectively. The 7r/6 bending refraction beams in the first two cases and the almost totally 
reflected beam in the third case are clearly observed. 

The self-collimation effects in square lattice PC or PhC have already been discussed 
by some researchers [TTH21"] . however, to our knowledge, this effect in triangular lattice 
structures has not yet been utilized as a means to design an acoustic lens for focusing 
sound. Here we show its applicability. Unlike most previously designed flat lens using 
circular EFC, we use a hexagonal EFC to achieve this goal. Usually a size-matched circular 
EFC having an effective refraction index close to —1 is preferred, because it has several 
advantages such as that for any incidence angle a refracted wave can always be excited, 
and all rays can be brought to the same focal spot. However, usually size-matched EFC 
can only be realized through sophisticated methods, such as using more than one kinds of 
host media [14]. Our aim in this paper is to design a simple PC lens without relying on 
those sophisticated methods, thus the hole-array structure is our choice. In order to reduce 
the total reflection power, the EFC size of the PC must be large enough and close to that 
of the host medium. This consideration leads us to choose ua/2nc t = 0.6 as the working 
frequency. If a slightly higher frequency EFC is chosen, the reflection gets higher, and the 
image intensity gets dimmer. The anisotropic refraction property of the PC looks like a 
disadvantage, but we will show that a flat lens of this PC can indeed work very well and it 
can brought the waves emitted by a point source to a image spot of subwavelength wide. 

Now we put a point source in front of the lens and study the imaging effect. Simulation 
results for 9-layer and 13-layer structures are shown in Fig. 3(a)-(c) and (d)-(f), respectively. 
From (a) to (c) or (d) to (f) the source-lens distance is increased from a to 6a. According 
to these simulation results, the shear waves inside the PC indeed propagate along specific 
directions (the ±7r/6) as expected. Besides, the thicker the lens, the farther the image 
locates, consistent with the negative refraction description. Furthermore, when the source 
moves towards the left, the image also moves towards the left, keeping the source- image 
distance almost unaltered. 

The intensities |w| 2 of the displacement field in the focal planes (located at the positions 
of the image peaks) for the 9-layer and 13-layer lenses are shown in Fig. 4(a) and (b), 
respectively. In each subplot three source-lens distances (d = 0.8a, 3.3a, and 5.8a) are 
examined, and the image widths w defined by their full-width-at-half-maxima (FWHM) are 
recorded. The results reveal that the lens can really make a subwavelength image of the 
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point source at the far field zone, although the image width (about 0.6A) is a little larger 
than that employing evanescent wave coupling mechanism [221 ES] ■ 

The position of focal point can be predicted by considering the refraction of acoustic rays. 
Although the PC in the chosen frequency cannot be represented by an isotropic effective 
medium, a formal usage of Snell's law (the consequence of the continuity of the tangential 
component of wave vector across the host-PC interface) can still help us to find the position 
of the focal point approximately The effective refractive index n = —0.8305 of the PC at 
reduced frequency k t a/2ix = 0.6 can be derived from the relation n = —c t \\t\/u} = — |k|/fc t , 
here k t is the radius of the host EFC, and |k| is the shortest radius (the line segment TP in 
Fig.l(b)) of the PC EFC. We sketch the ray-trace diagram in Fig. 5. The refraction angle are 
fixed at ±7r/6, as mentioned above. According to Snell's law, we have sin# = n sin(— 7r/6), 
which gives 9 = 0.4282 or 24.53° The distances x\ and 22 are related to d\ and d,2 through 
the relation X\jd\ = £2/^2 — tan 6/ tan(7r/6) = 0.7906. We thus have a constant source- 
image distant d± + c?2 + D = 2.265-D, consistent with the results obtained in the wave 
field simulations. This analysis seems to imply only one incident angle is allowable, which 
is not the case for waves emitted from a point source. The ray trace calculated in this 
way represents the route having an extremum phase aggregation, that is, waves interfere 
constructively along this specific route. If a route slightly deviated from this specific one were 
used, the phase propagation direction (the direction of k) and beam propagation direction 
(the direction of (S)) inside the PC would not be parallel to each other. This leads to 
partially destructive interference and reduce the strength of the ray. A simple glance of 
Fig. 3 makes these statements more concrete. The predicted positions for the two images 
(one inside and one outside of the lens) based on the ray formula and those obtained by 
wave simulations were compared in Fig. 6. A good agreement between these two methods 
is found, as can be easily checked. 

In conclusion, we have proposed a 2D phononic crystal flat lens embedded in solid matrix 
to achieve far-field subwavelength imaging of off-plane shear waves. The field patterns 
were computed by using the multiple scattering method, whereas the band structure of the 
phononic crystal were obtained by using the Dirichlet-to-Neumann map method. Image 
positions were derived by using ray-trace method, and the results are in good agreement 
with those from direct wave simulations. Our study has the following novelties. First, 
our PC lens focuses transverse shear wave, while most of previously proposed PC lenses 
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focus pressure wave. Second, although we choose aluminum as the host medium, the PC 
band structure is in fact universal when properly normalized, thus the same kind of lens 
can be fabricated using other solid matrix materials. Third, we demonstrated that the 
self-collimation effect of triangular lattice PC can be employed to focus elastic waves, thus 
circular shaped equal-frequency-contour is not necessary for achieving negative refraction or 
subwavelength imaging. Our study may provide some insights for designing new kinds of 
phononic devices. 
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FIG. 1: (a) The band structure of the triangular lattice hole-array PC. (b) The equal frequency 
contours (EFC) and the refraction diagram. The lattice constant is a, and the radius of the holes 
is ro = 0.4a. The two straight lines in (a) are the "sound lines" for the shear wave. The hexagonal 
blue curve marked 0.6 and the red circle in (b) represent the EFCs of the PC and of the host 
medium respectively, at reduced (dimensionless) frequency kta/2ir = u}a/2irct = 0.6. The two 
thick black arrows indicate the k vectors of the incident and refracted waves. The continuity of the 
wave phase along the host-PC interface implies that these two k vectors have the same component 
along the interface. The two pink arrows represent the energy flow (averaged Poynting vector) 
directions of the incident and refracted waves. The wave vector length |k| = TP is chosen for 
calculating the effective refraction index n. Note that the direction of the refracted ray (energy 
flow) is approximately fixed at ±7r/6. Besides, it is not possible to excite a propagating mode of 
the PC if the incident angle is larger than 7r/3. 

[35] Pi-Gang Luan, Kao-Der Chang, J. Nanophotonics 1, 013518 (2007). 
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FIG. 2: When a shear wave beam of (dimensionless) frequency 0.6 is incident upon the PC, it 
gets refracted with angle 7r/6 (if the incidence angle is smaller then 7r/3) or totally reflected (if the 
incidence angle is close to or larger then vr/3). The field patterns for incidence angles equal to 10°, 
40° and 60° are shown in (a), (b) and (c), respectively. 
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FIG. 3: The displacement field patterns under the illumination of a 2D point source. The center 
of the first layer of holes is located at x = a. The results for the 9-layer structure are shown in 
(a)-(c). The source is located at x = (a) 0, (b) —2.5a, and (c) — 5a, respectively. Similar simulation 
results for the 13- layer structure are shown in (d)-(f). 
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FIG. 4: The shear wave intensities in the focal plane for the (a) 9-layer and (b) 13-layer structures. 
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FIG. 5: Refraction of acoustics rays. The refraction angles of ±7r/6 are expected inside the lens. 
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FIG. 6: Comparisons between the predicted (the "+" sign) and the actual (the "o" sign) positions 
of the images. The vertical axis is the source-lens distance measured from the source location to the 
center of the first layer of the holes. The horizontal axis is the position of the image. The results 
for the 9-layer and 13-layer structures are shown in the upper and lower subplots, respectively. 
The two vertical lines represent the boundaries of the lens. 
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